Contents 



1 Introduction 



2 Symmetry as a Truncation Principle 






[2.1 Modern Approaches to the Shell Model Problem 




2.2 Symmetry-Dictated 'truncation 






2.3 Shell Model Effective Interactions . 


2.4 FDSM Effective Interactions 





3 Examples of Truncated Valence Spaces 



4 


Symmetry-Breaking in the FDSJVJ 


4.1 The Breaking of Dynamical Symmetries 




4.2 Symmetry Breaking by Spherical Single-Particle Energies 




4.3 Comments on the Size of Symmetry Breaking Terms 



5 Heritage-0 Calculations for Major Shell Truncation 

5.1 Hamiltonian 

5.2 Particle Number Distribution 

5.3 Masses of Heavy and Superheavy Nuclei: Z= 1 10 1 1 1| 





p. 3.1 Comparison of Mass Calculations for Very Heavy Elements 




p. 3. 2 Superheavy Elements 




p. 3. 3 Experimental Implications 


5.4 


The Xe-Ba Region 




p. 4.1 Energy Spectra 




p. 4. 2 Electromagnetic Transition Rates 




|5. 4. 3 Fermion Dynamical Symmetry versus Supersymmetry 




p. 4. 4 Tau Compression 


5.5 


SU(6)-Like Behavior in the Pt Region 




p. 5.1 Hamiltonian and Spectra 




p. 5. 2 Electromagnetic 'transitions 




p. 5. 3 F-Spin and Majorana Interactions 


5.6 


Effective Interaction for Rare Earth Nuclei 


5.7 


Universality of Normal and Exotic States 




p. 7.1 Hamiltonian and Operators 




p. 7. 2 Spectrum and Multipole Strengths 



6 Higher Heritage Configurations 



7 Conclusions 



Solution of the Nuclear Shell Model by 
Symmetry- Dictated Truncation* 



M. W. Guidry, 1 ' 2 D. H. Feng, 3 X.-W. Pan, 3 and C.-L. Wu, 4 

1 Department of Physics, University of Tennessee, 

Knoxville, TN 37996-1200 USA 

2 Physics Division, Oak Ridge National Laboratory 

Oak Ridge, TN 37831 USA 

^Department of Physics and Atmospheric Sciences, 
Drexel University, Philadelphia, PA 19104 USA 

4 Department of Physics, Chung Yuan Christian University, 
Chung- Li, Taiwan 32023 ROC 

—February 9, 2008— 
Abstract 

The dynamical symmetries of the Fermion Dynamical Symmetry Model are used as a principle 
of truncation for the spherical shell model. Utilizing the usual principle of energy-dictated 
truncation to select a valence space, and symmetry-dictated truncation to select a collective 
subspace of that valence space, we are able to reduce the full shell model space to one of manage- 
able dimensions with modern supercomputers, even for the heaviest nuclei. The resulting shell 
model then consists of diagonalizing an effective Hamiltonian within the restricted subspace. 
This theory is not confined to any symmetry limits, and represents a full solution of the original 
shell model if the appropriate effective interaction of the truncated space can be determined. 
As a first step in constructing that interaction, we present an empirical determination of its 
matrix elements for the collective subspace with no broken pairs in a representative set of nuclei 
with 130 < A < 250. We demonstrate that this effective interaction can be parameterized in 
terms of a few quantities varying slowly with particle number, and is capable of describing a 
* Review article prepared for the Journal of Physics G. 
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broad range of low-energy observables for these nuclei. Finally we give a brief discussion of 
extending these methods to include a single broken collective pair. 
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1 Introduction 

The shell model is fundamental to the understanding of nuclear structure, but it cannot be 
used for practical calculations in systems with many valence particles outside closed shells. 
There are two basic reasons for this. (1) The matrix dimensionalities are too large. (2) There 
are too many effective interaction parameters. The first problem is well known; the second is 
as important, but is less appreciated. We may illustrate this second problem by noting that 
if one views the effective interaction for a major shell of neutrons and of protons as specified 
by a set of matrix elements to be determined by fits to existing data, for the heaviest nuclei 
we must determine a minimum of about 2500 parameters (matrix elements of allowed one 
and two-body interactions) from some combination of theory and data. By comparison, the 
same approach in the sd shell requires almost two orders of magnitude fewer parameters. 

2 Symmetry as a Truncation Principle 

This proliferation of parameters is not a failure of the theory. The large number of parameters 
just reflects the number of independent one and two-body scatterings that are possible within 
a large valence space populated by significant numbers of nucleons. Of course in all practical 
calculations, the number of parameters has been reduced by some means. The most drastic 
prescriptions are the mean field ones that replace these parameters by an average potential, 
and an additional set of prescriptions (often in the form of semi-empirical recipes) that 
prescribe how to calculate observable quantities, and how to correct at some level (again 
often semiempirical) for the correlations neglected in such an approximation. 

The applicability of mean field theories in a variety of applications is well established, 
but the nuclear many-body problem cannot be reduced to mean field terms, and the semiem- 
pirical nature of mean field successes argues for a more microscopic understanding. On the 
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other hand, the success of mean field methods, and the relative simplicity of much of low-lying 
nuclear structure, suggest that most of the possible contributions to the effective interaction 
either are negligible, or more likely, they enter (low energy) nuclear structure only in certain 
restricted combinations. 

2.1 Modern Approaches to the Shell Model Problem 

Thus, the solution of the general shell model problem for heavy nuclei requires a resolution 
of the matrix dimensionality problem and a consistent method to select a highly restricted 
subset of the effective interaction parameters as the ones relevant for low energy structure. 
Ideally, one would like an approach that accomplishes both of these tasks in a self-consistent 
manner. Traditionally one begins with some form of energy-dictated truncation. This helps 
with the matrix dimensionality problem by restricting the size of the shell model space, and at 
the same time limits the number of parameters that must be determined. However, energy- 
dictated truncations have had limited success for nuclei far removed from closed shells: 
the influence of any single high-lying configuration may be negligible, but the aggregate 
contribution of many such configurations may not be. Thus, energy-dictated truncations 
alone are unlikely to allow us to solve the general shell model problem for heavy nuclei. 

There are four modern approaches to this problem of trying to extend the shell model 
to the description of heavy nuclei far removed from closed shells. 

1. Improved algorithms and computers for traditional shell models (e. g., Refs. |1|, |2], |3|). 

2. Path integral solutions of the shell model using Monte Carlo algorithms on fast super- 
computers (e. g., Refs. ||, 

3. Truncations of the shell model space based on guidance from mean-field geometrical 
models with approximate pairing (e. g., Refs. || [7|]). 

4. Symmetry-dictated truncations of the shell model space (e. g., Refs. || |9], [10], |TT| ) . 
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Approaches (1) and (2) have had some success with the matrix dimensionality problem, 
either by attacking it more efficiently, or by avoiding it altogether in the path integral 
approach. However, neither offers an intrinsic solution to the burgeoning number of effective 
interaction parameters in the heavy nuclei. One still must supplement these approaches with 
a prescription for selecting the components of the interaction that are to be emphasized. We 
will not address methods (1) or (2) in this review, and refer the reader to the literature cited 
above on these subjects 0, |j, |], 

On the other hand, approaches (3) and (4) help with both aspects of the problem: 
the truncations implied by these prescriptions reduce the matrix dimensionality, and the 
nature of the truncation suggests a prescription for choosing the form of the most important 
interaction terms in the resulting trunctated space. An excellent example of approach (3) 
is the Projected Shell Model 0, which truncates the space using a deformed mean field 
plus BCS pairing prescription, and then diagonalizes a Hamiltonian composed of low-order 
pairing and multipole terms in the truncated basis. The approaches in category (3) have 
had considerable success in practical calculations, but we will not deal with them in this 
review, except to remark that they have many deep similarities with the symmetry-based 
approaches that we shall discuss. 

We focus the remainder of our discussion on symmetry-dictated truncation and its 
integrated solution to the shell model problem: (1) the symmetries dictate a severe trun- 
cation of the shell model space; (2) the requirement that the dominant interactions respect 
these symmetries provides a methodology for emphasizing a limited subset of interactions. 
Whether such an integrated approach is useful depends on the ability of the resulting theory 
to describe a broad range of nuclear structure data with a restricted set of parameters that 
vary in a well-conditioned manner with particle number. 
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2.2 Symmetry- Dictated Truncation 

Symmetry-dictated truncation is illustrated for a simple model assuming an SU2 symmetry 
in Fig. [I], and for a shell model dynamical symmetry in Fig. |^. Stated somewhat loosely, an 
energy-dictated truncation truncates the space "spherically" in the space of symmetry gen- 
erators (for example, the three angular momentum components in Fig. [I]), but a symmetry- 
dictated truncation reduces the space further by selecting a particular "direction" (or set 
of directions) in the space of symmetry generators. Representative directions are indicated 
schematically in Fig. [I] and Fig. ||] by the heavy arrows. Such an approach is closely associ- 
ated with spontaneous symmetry breaking and phase transitions because of the selection of 
preferred directions in the space. 

2.3 Shell Model Effective Interactions 



Shell model effective interactions are highly dependent on the truncation scheme |13], |14 , 



15 1 . A realistic initial approach to obtaining an effective interaction is to treat its matrix 



elements as parameters to be determined from data[16|. The prototype of this method is the 



determination of the Brown and Wildenthal matrix elements for light nuclei [|T7j |18| , but a 
literal transcription of the empirical sd shell approach to the actinide nuclei would require 
approximately 2500 matrix elements of the effective interaction to be determined from the 
data. 

On the other hand, in the actinide region the most general Sp$ x Sp$ Hamiltonian of the 
Fermion Dynamical Symmetry Model || requires (at most) about 100 effective interaction 
parameters to be determined, which is comparable to the number required in the full shell 
model for the sd shell ||17|| . In reality, the symmetry-dictated truncation of the FDSM will 
limit the number of relevant parameters even more severely than these formal estimates would 
suggest: the results to be presented here indicate that most low-lying nuclear properties are 
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determined by only of order 10 symmetry-selected parameters for the entire shell, with (at 
most) a weak dependence of these parameters on particle number. 

2.4 FDSM Effective Interactions 

Thus, the FDSM symmetry-dictated truncation provides a methodology with the potential 
to enable shell model calculations for all nuclei. To establish its validity, it is necessary 
to demonstrate through a detailed set of calculations that the required effective interaction 
(1) is easily determined, and (2) produces a quantitative description of observables in heavy 
nuclei. This will be an iterative process requiring a long series of systematic calculations, 
but we offer the following initial speculations concerning the FDSM effective interactions and 
the process of determining it: 

1. The symmetry-dictated truncation of the FDSM is qualitatively different from energy- 
dictated truncations; it could lead to effective interactions in the highly truncated 
spaces that are markedly different in the components that are emphasized relative 
to the appropriate interactions in the full shell model space. Unlike energy-dictated 
truncations, a symmetry-dictated truncation selects particular "directions" in the shell 
model space (see Fig. Q). 

2. A complicated particle number dependence for the effective interaction would make the 
theory difficult to use in a consistent microscopic way. Therefore, we require that the 
effective interaction have a smooth and/or weak particle number dependence within 
a shell. Since filling different valence shells defines different phases of the theory, it 
is possible (but not required) that the parameters could have discontinuities between 
shells. Likewise, it is possible (but not required) that there may be discontinuities 
between parameter sets describing regions having different dynamical symmetries, since 
these too define different phases of the theory. 
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3. A reasonable first step will be to determine the effective interaction in the SD collective 
subspace of the FDSM using symmetry-limit calculations. We may then investigate 
deviations from the symmetry limits within the SD space, and may expect that this 
leads to small and smooth changes of the effective interaction parameters from the 
symmetry- limit values. 

4. Finally, we may expand the space to include broken collective pairs. This will allow the 
determination of additional effective interaction parameters associated with high-spin 
physics that played small roles in the low-lying states. In addition, we may expect a 
systematic renormalization of the parameters already determined in the SD subspace 
by virtue of enlarging the space. This renormalization under space enlargement must 
be well behaved if the symmetry-dictated method of truncation is to be of practical 
utility. 

Thus an important first step in this program is to take dynamical symmetry limits as a 
starting point and fit to the data using an FDSM Hamiltonian with symmetry breaking 
terms included in lowest order perturbation around the symmetry limits. Applications of the 
FDSM at this level have been discussed extensively in ||. At the next level of improvement, 
we may assume the most general highest symmetries and a corresponding Hamiltonian, but 
not restrict the calculation to the dynamical symmetry limits or to perturbation around the 
symmetry limits. We may hope to obtain by these approximations a qualitative agreement 
with experiment; this effective interaction can then be refined in subsequent calculations that 
include more realistic configuration mixing and enlarge the space to include broken pairs. 

In this review we shall emphasize two aspects of this extension of the FDSM beyond the 
dynamical symmetry limits. The first is the initiation of a systematic analysis of symmetry- 
breaking terms in the realistic shell model Hamiltonian relative to the FDSM symmetry 
limit Hamiltonians. The second is a series of numerical calculations that take the theory 
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beyond the symmetry limits and begin to establish a set of effective interaction parameters 
for realistic symmetry-truncated shell model calculations. 

3 Examples of Truncated Valence Spaces 

The essence of the FDSM approach is to impose a symmetry-dictated truncation within a 
particular valence space, with the choice of valence space often dictated by energy consider- 
ations. As has been discussed in some detail in Ref. ||, there is considerable flexibility in 
the choice of the space in which to implement such a truncation. The space truncation is 
instituted through the generalized Ginocchio [|B| coupling scheme illustrated on the far right 
side of Fig. [3] f|. The most common implementation has been within a single major shell, 
as illustrated by the coupling labeled "Normal Deformation" in Fig. || However, an FDSM 
model of superdeformation has also been proposed that uses the alternative coupling 
scheme labeled "Super Deformation" in Fig. |3|. This example makes it clear that the choice 
of valence space is not the central issue, but a clever choice of the valence space may enhance 
the ability of the symmetry limits of the theory to describe data. 

More generally, it has been suggested that the FDSM can be formulated in a variety 
of valence spaces, with different choices associated with what would be termed shape co- 
existence in a mean-field theory pi| . In this review we shall concentrate on applications 



of symmetry-dictated truncation to major-shell valence spaces, but the basic ideas should 
be applicable to a broad range of valence spaces if one can implement a symmetry-dictated 
truncation within them. In particular, we expect that a similar analysis is possible for 
superdeformed states. 
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4 Symmetry-Breaking in the FDSM 

The symmetry limits of the FDSM, and perturbation theory about those symmetry limits, 
describe a broad range of nuclear structure observations ||. It is of interest to examine 
systematically the excursions from the symmetry limits of the theory, in order to test its 
suitability as a systematic truncation procedure for quantitative shell model calculations in 
heavy nuclei. 

4.1 The Breaking of Dynamical Symmetries 

Principles of dynamical symmetry and dynamical symmetry breaking are summarized in 
Fig. f|. A dynamical symmetry results when a Hamiltonian can be expressed in terms of 
invariants for the highest group and subgroups of a group chain, as we illustrate for the 
Sp$ D SU3 D SO3 dynamical symmetry of the FDSM. Symmetry-breaking then corresponds 
to the presence of Hamiltonian terms that cannot be expressed in this manner. Figure ^ 
also distinguishes two approaches to dynamical symmetry breaking. The first is empirical: 
one defines the symmetry breaking to be the difference between the observation and the 
symmetry- limit calculation. This provides a definition of symmetry breaking, but is not 
predictive. More useful are the microscopic approaches to symmetry breaking, where the 
symmetry- limit Hamiltonian H sym is presumed to derive from a more fundamental Hamilto- 
nian -ff m i cr o with a known form. Therefore, one can use theoretical and physical arguments 
to identify likely symmetry breaking terms if s b, and can make predictive estimates for the 
expected magnitude of symmetry breaking. The situation addressed here falls in this second 
category: the FDSM symmetry-limit Hamiltonian represents an approximation to the full 
shell model Hamiltonian that omits terms breaking a particular dynamical symmetry. Since 
the forms of the two Hamiltonians are known, we may use physical arguments to identify 
the terms of the full Hamiltonian that break the symmetry and are likely to be important. 
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We now illustrate using Sp G D SU 3 symmetry breaking by spherical single-particle energies 



in the FDSM 22 



4.2 Symmetry Breaking by Spherical Single-Particle Energies 

Orbitals exhibiting an Sp(6) D SU(3) dynamical symmetry have k = 1 (see Fig. |||). A 
Hamiltonian with an Sp(6) D SU (3) dynamical symmetry requires degeneracy in the single- 
particle energy terms corresponding to the same value of i. For normal deformation in 
heavier nuclei, and for super deformation, there are multiple values of i within a valence 
shell and the symmetry-limit Hamiltonian will exhibit a higher level of degeneracy than the 
realistic spherical single-particle shell model spectrum (Fig. [3]). The difference between the 
symmetry- limit and realistic spectra of Fig. [3] are shown in Fig. [|, which also summarizes 
several issues that are relevant to the present discussion. These concern, not simply the size 
of the single-particle splitting, but its size relative to the correlation energy of the system 
and how much of that splitting breaks the relevant symmetry |50, |8|, 2T| . 



The states of the FDSM are classified according to a total heritage quantum number u 
that measures the number of particles not coupled to coherent S and D pairs [fU ]. Low-spin 



states of even-even nuclei are dominantly u = configurations. The mixing matrix elements 
associated with the splitting of the single-particle energies may be expressed as 



X'n'u' 



n 



(rr)0 



Xpu 



(X'f/u' 


[nf r)0 , Csu(3)} 


Xfiuj 


c(\fi) - c(Ay) 





where u and u' are the heritage quantum numbers, C[Xp) is the usual eigenvalue of the 
quadratic SU3 Casimir operator Csu{S) evaluated in an SU3 representation (A,/i), and the 
single-particle energy expressed in terms of the standard FDSM k-i basis is: 



EV^ (rr)0 r 
j r,i 



(2) 
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where the operator n'f^ is given by 



n 



(rr)0 



2a 



h lhi 



(rr)0 



(3) 



and 



k i j 
k i j 
r r 



(4) 



with the square bracket denoting a normalized 9-j coefficient. The k-i basis b ki has been 



defined in 



Qj and fij are the pair degeneracies for the j shell and the shells 



associated with pseudospin i, respectively [flj = j + |, and = (2k + l)(2i + l)/2]. Ref. 



22| , demonstrates that the only operator capable of mixing an SU3 irrep in the u = bands 



with other irreps is n,- 11 ' )0 . Estimates for the upper limit on the symmetry breaking 
associated with this term are summarized in Fig. |6] for typical heavy, normally-deformed 
nuclei. A similar analysis for superdeformation suggests that the single-particle symmetry 
breaking associated with the FDSM model of superdeformation is even smaller than that 
exhibited in this example for normal deformation. Thus, we expect that spherical single- 
particle symmetry breaking in both normal and superdeformed nuclei will be perturbative 
in size for the Sp^ Z> SU3 symmetry of the FDSM, and the full inclusion of such terms will 
be unlikely to invalidate previous symmetry-limit results. 

4.3 Comments on the Size of Symmetry Breaking Terms 

The stability of the SU3 dynamical symmetry for the FDSM derives from the particular 



structure of the Ginocchio S-D pairs |T9|, |1J, ||. The single-particle terms can break the 
SU3 symmetry only through an indirect Pauli effect associated with the embedding of SU3 
in the higher symmetry Sp 6 . The only operator in the single-particle energy terms that 



can accomplish this is n. 



(11)0 



In other fermion theories, such as the Elliott model 



or 
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the pseudo-S , [/ 3 model [£5|, [26J, the SU 3 symmetry is embedded in a much larger group and 
there are many generators that could break the symmetry directly 

This illustrates an important concept concerning symmetry breaking: the size of the 
symmetry-breaking terms for a particular dynamical symmetry may differ considerably from 
that expected based on experience with mean field theories, or even theories based on a 
symmetry that is formally related but implemented in a different physical manner. In this 
case we see that the single-particle terms have little influence on the SU3 symmetry because 
of the unique properties of the pairs that generate the symmetry. On the other hand, in a 
deformed mean field theory one generally finds that the quantitative results are much more 
sensitive to the details of the initial spherical single-particle spectrum. 

5 Heritage-0 Calculations for Major Shell Truncation 

Let us now examine some FDSM calculations. The majority of these will employ the code 



FDUO of Wu and Vallieres p7| , [28] that diagonalizes the most general interaction consistent 
with the highest shell symmetry of the valence space. Thus, its solutions may be considered 
to be linear combinations of the dynamical symmetries allowed in a valence space. It is valid 
only for heritage-0 subspaces; thus, it is most applicable for states in even-even nuclei below 
angular momentum 10, where the breaking of pairs is not too important. In addition, this 
code assumes that N\, the number of pairs in the normal parity orbitals of the valence shell 
(for neutrons or protons), is conserved. Thus, it does not incorporate directly the effects 
of scattering particles between the normal and abnormal parity orbitals (such effects may 
still be included indirectly in the effective interaction [^, 0). Because this scattering is 
expected to be more important in the Pauli forbidden region lying between 1/3 and 2/3 
filling of the normal parity orbitals in Sp^ shells H , the validity of FDUO in that region 
is also questionable. Therefore, in this review we shall confine attention to states of low 
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angular momentum for nuclei that do not lie in the Pauli forbidden region of Spe shells. For 
SO$ shells, there is no such restriction and we are free to apply FDUO calculations to all 
nuclei in the shell. 

5.1 Hamiltonian 

The most general Hamiltonian consistent with an SOs or Sp^ highest shell symmetry may 



be expressed as [27, 



+B^P^-P 2 U + B r p 3- p 3 + B i p i- p i + B\Pl-Pl + B™P»-P? (5) 

where 5* denotes monopole pair operators, D denotes quadrupole pair operators, and P r 
denotes multipole operators of order r, with all quantities constructed in the k-i truncation 
scheme illustrated in Fig. |3|. Not all of these terms contribute for a particular highest 
symmetry. One finds that for spectra (which depend only on energy differences) there are at 
most 11 parameters for SO$ x SO$, 8 for Sp^ x Sp$, and 9 for SO^ x Sp$. In the numerical 
calculations that follow, we will typically use a simplified version of this Hamiltonian that 
retains only pairing and quadrupole terms, thereby reducing the number of free parameters 
to 5 or less for most cases. 

5.2 Particle Number Distribution 

In the calculations to be presented in this review, the number of pairs N\ in the normal- 
parity levels is treated as a good quantum number and is estimated from the semi-empirical 
formula determined globally from the ground state spin of the odd-mass nuclei ||: 

N for iV < 1.5 

Nx = { 0.75 + 0.5N for 1.5 < N < 2Q Q + 1.5 (6) 
N -Vl for N > 2fi + 1-5 
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where N is the number of valence pairs and fio is the pair degeneracy of the abnormal-parity 
level. The above formulas are suitable for both protons and neutrons, and are found to be in 
excellent agreement with a similar distribution obtained empirically from the Nilsson level 
scheme with a measured deformation parameter (3. 

5.3 Masses of Heavy and Superheavy Nuclei: Z=110— 111 



In recent years, evidence has accumulated |32|, [33], for extension of the known elements to 
larger proton number, culminating in the discovery of several isotopes having proton number 



Z = 110 — 111 ||35|| . The usual view is that these isotopes represent the tail of the "normal 
element" distribution, and that the predicted superheavy elements represent a qualitatively 
different set of nuclides that remain undiscovered at a predicted proton number of about 
114 and a neutron number approximately 20 units heavier than for the recently-discovered 
isotopes of elements 110 and 111. As we now discuss, truncated shell model calculations 
suggest a different interpretation of these results that could have important experimental 
implications |!6], j3Tf . 

5.3.1 Comparison of Mass Calculations for Very Heavy Elements 

In Table [I], we compare the observed mass excesses (first row, labeled "Exp" ) for isotopes of 
the heaviest elements with calculations for these masses. A perusal of these results suggests 
that the FDSM mass calculations described in Ref. [|36| give a global description of the 



heaviest elements that is arguably the best now available. This is particularly noteworthy 
because the FDSM calculations were not tuned specifically to these heaviest elements; the 
parameters employed in Ref. [|36| were determined by a global fit to all available masses above 
Z = 82 and N = 126. We also emphasize that the accurate description of the Z = 110 — 111 
isotopes represents a pure prediction of the FDSM theory, published before their discovery, 
and without adjustments to the theory based on specific properties of the Z = 108 — 109 
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elements that were known at the time of the calculations. 
5.3.2 Superheavy Elements 



An island of superheavy elements was found in the FDSM mass calculations ||36|| , but this 
island was shifted to considerably lower neutron number than is predicted in classical cal- 
culations of the stability of superheavy elements, was found to be more stable than in such 
calculations, and was found to correspond to nuclei having near-spherical shapes. The shell 
correction associated with this minimum is illustrated in Fig. [7] for the case of a spherical 
Woods-Saxon single-particle spectrum. The location of the new predicted maximal shell sta- 
bilization, the traditional location of the island of superheavy stability, and the new isotopes 
of elements 110-111, are also shown in Fig. 0. These properties of the FDSM heavy-element 
solutions are consequences of two general physical principles. The first is a monopole- 
monopole interaction that is expected to be present on general shell-model grounds, and 
found from empirical mass fits to become increasingly repulsive in moving away from the 
stable nuclei. The second is that the structure of these heaviest elements is expected to 
involve a competition between an SU2 dynamical symmetry favoring spherical symmetry 
and an SU3 dynamical symmetry favoring axially symmetric deformation [||]. 

5.3.3 Experimental Implications 



Since Ref. |36| was a general treatment of nuclear masses that was not tuned to predict the 
properties of these heavy new isotopes, we are emboldened to draw some general conclusions 
from this successful prediction that differ from those of more traditional approaches to the 
stability of the heavy elements 

1. The observation of isotopes with proton number Z — 110 — 111 having the predicted 
properties may indicate a shift of the superheavy island of stability to considerably 
lower neutron numbers, as predicted in Ref. |36] . 
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2. These new isotopes are expected to be near-spherical or deformation-soft structures, 
in contrast to the heavy elements with Z < 106, which are found in the FDSM (and 
other) calculations to have well-deformed ground states, and to most other descriptions 
of the new Z — 110 — 111 elements, which find them to have deformed ground states. 

3. The successful prediction of the masses of these new isotopes is further evidence sup- 



porting the importance of the monopole-monopole interaction introduced in |36|] for 
the description of nuclear masses. 

4. No shell stabilization minimum is expected at the traditional location of the superheavy 
elements because of the repulsive monopole-monopole interaction. Thus, we predict 
that the traditional island of superheavy elements does not exist, and that the r-process 
element production path may run closer to the stability valley than has generally been 



assumed for the heaviest elements 38 



The present hypotheses have some consequences that are testable, though the experiments 
are difficult. (1) There should exist nuclides of elements Z = 112-114 that are approximately 
as stable as those of elements 110-111 (see Fig. [j]). (2) Beyond Z ~ 116 and N ~ 170, the 
heavy nuclei should rapidly become less stable, and the region of traditional superheavy 
nuclei should be quite unstable (see Fig. [I]). (3) The nuclides in this new region of super- 
heavy nuclei (Z ~ 110 — 116 and N ~ 160 — 170) are expected to be either spherical or 
very deformation soft, with attendant consequences for properties such as the alpha decay 
systematics. (4) The implied shift of the r-process path closer to the stability valley should 
also have observable consequences, but this may require an improved understanding of the 
astrophysical environment for the r-process. 
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5.4 The Xe-Ba Region 

The FDSM is defined by a fermionic Lie algebra, has symmetry limits analogous to all the 
IBM limits, and takes the Pauli principle into account fT9| , R{|. Furthermore, the states for 
even and odd systems in the FDSM belong to vector and spinor representations, respectively, 
of SOg or Spe, thus allowing a description of even-even and even-odd nuclei without the 
necessity of additional degrees of freedom. Therefore, symmetry-truncated shell model calcu- 
lations based on the FDSM dynamical symmetries for even and odd nuclei offer a microscopic 
alternative to IBM and nuclear supersymmetry (NSUSY) as means to unify the structure of 
even and odd systems. 

Nuclei in the Xe-Ba region have neutrons and protons in the 6th shell, with FDSM 
pseudo-orbital angular momentum k = 2 and pseudo-spin i — | for the normal-parity levels. 
Thus they are expected to possess SOg x SO$ symmetry, which contains an SOg +u subgroup, 
and to have analytic solutions for 5*05 x SU2, SOq, and SO7 dynamical symmetries. For 
these reasons, the Xe-Ba region is an excellent one in which to investigate a unified FDSM 
description of even-even and even-odd nuclei. 

5.4.1 Energy Spectra 

The wavefunctions for both even and odd nuclei are given by the following group chain 

( S0\ D SO\ D SOl)x SO k 5 D SO l 5 +k D S0% +1 (7) 
[hhhk] [0-10-20-3] [tit 2 ] [t] [UlU) 2 ] J 

where [hhhh], [010203], and [tit 2 ] are the Cartan-Weyl labels for the groups S0$, SOq, and 
SO5, respectively, r = 0(1) for even (odd) nuclei, and k and i indicate pseudo-orbital and 
pseudo-spin parts of the groups, respectively. The energies for even systems are 

E even = + g 6 cr(a + 4) + g' 5 r(r + 3) + g'jJ(J + 1), (8) 
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and for odd systems, 

^odd = E (o) + g^fa + 4) + ^ + 2 ) + (a 3 ) 2 } + 9J J(J + 1) 

+ (#5 " ftObifa + 3) + r 2 (r 2 + 1)] + g 5 [wi(wi + 3) + cj 2 (cu 2 + 1)]. (9) 
More details on the construction of these formulas, and the reduction rules for even and odd 



systems, are given in Ref. [40 



The low-energy spectra for 120_132 Xe isotopes predicted by Eq. (Q) are compared with 
data in Fig. || and Fig. |9], and the parameters used in the calculations are given in Table [| 
The experimental spectra indicate that the SO3 parameter g\ is not sensitive to the neutron 
number in fitting the spectra of a chain of isotopes (including both even-even and even-odd 
nuclei). Therefore, in fitting the even-even nuclear spectra we fix the g\ parameter to be 
11.9 keV and the adjustable parameters were taken to be g§ and g' 5 . 

Using Table and Eq. (S), the spectra of the neighboring even-odd Xe isotopes can be 
calculated. We present the calculated and experimental results for 127 ~ 133 Xe in Figs. [TTJ [Til 
with the parameters given in Table [| Similar results for Ba isotopes may be found in Ref. 
40|| . Apart from the constant term, the formula for E even contains three parameters, and 



that for E odd contains two parameters beyond the three parameters that are determined 
by fitting the spectra of neighboring even-even nuclei. With two extra parameters, Eq. (§) 
reproduces the spectral patterns for the nuclei 127 ~ 133 Xe and 131_135 Ba (not shown here; see 
Ref. [[40]) with about 15 levels each. 



5.4.2 Electromagnetic Transition Rates 



In Ref. [40[, formulas are presented for the reduced quadrupole transition strengths B(E2) 
in both even and odd systems. Table [5] summarizes some B(E2) values calculated for even- 
even Xe isotopes. Similar results for Ba isotopes may be found in Ref. El. These are seen 



to agree rather well with the observations, and are indicative of a good SOq symmetry in 
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this region. In Table || both experimental and theoretical B(E2) values for the even-odd 
nuclei of 129 Xe and 131 Xe are given, and compared with the results of NSUSY calculations. 
The agreement with data is comparable in the two cases. 

5.4.3 Fermion Dynamical Symmetry versus Supersymmetry 



In Ref. it is shown that when the u = 1 fermion SOq state for the FDSM core is replaced 
by a boson state and the Pauli factor is ignored, the FDSM wavefunction goes over to an 
NSUSY wavefunction. Thus, NSUSY can be obtained as an approximation to the FDSM for 
odd-mass SOq nuclei. The full FDSM without this approximation implies corrections to the 
NSUSY picture. 

In the £7(6/4) NSUSY, the ground-state representation is [||] and the ground-state spin 
is always |. However, both | and | are observed as ground-state spins in this region. The 
second SO5 Casimir operator in Eq. (§) allows this possibility: for alternative signs of the 
parameter g 5 , the ground state spin can take the values ~ or | . What is more, by allowing 
#5 to change smoothly from positive to negative, the systematic shift of the ground band for 
the Xe (and Ba) isotopes can be reproduced, as shown for Xe in Fig. [12|. 

The general conclusion of Ref. ||(J is that the FDSM symmetry-dictated solution in the 
mass-130 region provides a unified description of even and odd nuclei that is comparable to, 
or better than, IBM and NSUSY approaches, but with a firmer microscopic basis. As one 
example of this microscopic basis, the group chain corresponding to odd-mass S0 6 nuclei 
in the FDSM is very similar to the corresponding NSUSY group chain, but pseudo-orbital 
angular momentum 2 is an ad hoc introduction in the NSUSY, while in the FDSM it is a 
natural result of the k-i reclassification for the shell model single-particle states in the sixth 
major shell. 
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5.4.4 Tau Compression 

Pure 5*06 spectra give reasonable agreement with data for r < 3 states in 120_126 Xe (see Figs. 
^|-|9|). However, the experimental energies for the higher r values in 128_132 Xe are lower than 
predicted, with the discrepancies increasing with increasing r. Energy levels within the same 
SO5 irrep (same r) follow the J(J + 1) rule rather well, so this discrepancy cannot be caused 
by the usual stretching effect. In Ref. [41|, it was shown that this is an SO5 r-compression 



effect whose driving force is the reduction of pairing correlation with increasing r. Allowing 
for g s = G — G 2 7^ 0, thereby deviating from the S0 6 limit, and treating the g s S^ ■ S term 
as a perturbation leads to the following energy formula, 

ELen = 4° + 9^ + 4) + A't(t + 3) - B'[t(t + 3)] 2 + g' r J(J + 1). (10) 
Fig. 13 shows the spectrum for 126 Xe predicted by Eqns. and (10), respectively. Inclusion 



of the SO5 stretching effect improves the agreement with data significantly. Details and 



more examples may be found in Ref. |41 . 



5.5 SO(6)-Like Behavior in the Pt Region 

Numerical studies using the FDSM have shown that there is an effective SO(6) fermion dy- 
namical symmetry in the 196 Pt mass region ||2 , [|3|] . More recently, a systematic description 



of nuclear structure has been given in this region using symmetry-truncated shell model 



methods [43[ . In this section, we review these calculations for the platinum isotopes and the 
emergence of effective 50(6)-like behavior from numerical calculations near 196 Pt, where the 
FDSM has no formal SOq symmetry. 

5.5.1 Hamiltonian and Spectra 

The five-parameter Hamiltonian is taken to be 

H = G'^SlS^ + G' 0u StS u + B^F* ■ P 2 n + B' 2v P 2 v • Pj + B 2nu P* • P 2 . (11) 
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In the SO (8) and Sp(Q) algebras, the quadrupole pairing interaction in this Hamiltonian is 
not independent and can be accounted for by redefining the parameters G' 0a = G Qa — G 2(7 
and B' ra = B ra — G 2(7 (c = ^, v). The spectra of five even platinum isotopes 19 °- 198 pt (there 
are approximately 13 levels for each nucleus) are found to be satisfactorily described by sets 
of five parameters. In Figs. |14] and [15] the calculated spectra for 190 ~ 198 pt are compared with 
the data; the corresponding parameters are summarized in Table [7]. 

Fig. [15] indicates that most low-lying + and 2 + states are well reproduced by the FDSM 
calculations, tolerable agreement is obtained for the 4 + states, and there is a persistent 
deviation of the calculated 3 + states that increases with mass. The basic theoretical trends 
are the narrowing of the energy gap between the 3^ and 4^ levels, and a widening of the 
gap between the 3+ and 2\ levels. This level pattern is characteristic of the SOq dynamical 
symmetry and appears naturally in the heavier Pt isotopes. In Fig. [14], we have presented 
a separate spectrum for 196 Pt, in order to show the experimental E2 branching ratios as 
well as the corresponding FDSM and IBM-2 predictions. The overall FDSM results are in 
quantitative agreement with the 0(6) level pattern and the SOq E2 selection rules in the 
FDSM description of 196 Pt. This agreement confirms that the FDSM does not have an explicit 
mathematical 5*0(6) dynamical symmetry for this case, but it has a remarkably accurate 
practical one j|2] . 



5.5.2 Electromagnetic Transitions 

The FDSM wavefunctions obtained through the diagonalization of the Hamiltonian can be 
used to compute multipole transition strengths. Expressions for quadrupole moments and 
B(E2) values, dipole moments and g-factors, and isomer and isotope shifts may be found 
in Ref. |43|. In Table H we have compiled the experimental B(E2) values and predictions 



from the FDSM and various versions of the IBM for 196 Pt. Additional comparisons for the 



22 



M. W. Guidry, D. H. Feng, X.-W. Pan, and C.-L. Wu 



other isotopes of Pt may be found in Ref. [ 43| . Also shown in the Pt table are SOq 



limit results. (IBM-l(g) means IBM-1 with one g-boson explicitly included and the effective 
charge individually fitted from nucleus to nucleus; in the IBM-2 calculation the effective 
charge remains constant for all the platinums, just as for the FDSM calculations.) Table |9] 
compares experimental g-factors with FDSM and IBM calculations, and calculated isomer 



and isotope shifts are compared with data in Table |10 
5.5.3 F-Spin and Majorana Interactions 

In boson-model calculations, F-spin is used to distinguish low-lying symmetric and mixed 
symmetry states [03], and a "Majorana interaction" is typically introduced to guarantee that 



mixed symmetry states do not appear at unacceptably low energies. For fermions there is no 
closed F-spin algebra and no explicit Majorana interaction in the Hamiltonian: for a fermion 
many-body system, the relative positions of the symmetric and mixed-symmetry states are 
governed by the basic pairing and multipole interaction strengths. For the 196 Pt case that 
we discuss, the mixed-symmetry 1 + state is at an excitation energy of approximately 2 MeV, 
and is correctly predicted without introducing the analog of a Majorana term. Thus, the 
effective fermion SO(Q) symmetry produces a quantitative agreement with the spectrum 
using a simpler Hamiltonian than is required in the IBM-2. 

5.6 Effective Interaction for Rare Earth Nuclei 

The parameters used in reproducing the 50(6)-like spectrum and the branching ratios of 
the platinum isotopes do not represent a peculiar set chosen solely for this purpose. Similar 
calculations in the remainder of the rare earth nuclei indicate that the Pt parameters are re- 
lated to the effective interactions of the entire shell. For example, if the same five parameters 
are used in the FDUO calculation but the appropriate valence particle numbers for heavy Gd 
isotopes are chosen, the spectrum and transition rates become SU (3)-like (axially-symmetric 
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rotors). With only small adjustments in these parameters (for example, less than 5% for the 
n-p quadrupole coupling, which is the dominant parameter) the spectrum and transition 
rates for the low-lying states in the heavier Gd isotopes can be quantitatively reproduced. 
Likewise, a similar set of parameters leads to vibrational (SU (2)-like) spectra for valence 
particle numbers near the beginning or the end of the shell. 

These results imply that a fixed Hamiltonian with parameters varying slowly with 
particle number can produce spectra that evolve from 577(2) to SU(3), from SU(3) to 
50(6), and from £0(6) to SU(2), as the proton or neutron number changes. This SU(2)- 
SO(Q)-SU(3) triangle is similar to the S77(5)-0(6)-S77(3) triangle relation [g|] of the IBM. 
However, unlike the IBM, the triangle relationships in the FDSM are shell dependent. For the 
50(8) x 50(8) shells, there is no SU(3) limit and the 50(6) symmetry is the most collective 
limit (corresponding to 7-unstable rotations). In this case the triangle relationship is replaced 
by 5£/(2)-50(7)-50(6). We have speculated J|2| that whenever there is a vibrational 



(SU(2)— like) symmetry limit on one corner of the triangle and an SU(3) — like rotational 
limit on a second corner of the triangle, an 50 (6)-like effective transitional symmetry may 
occur, even though there is formally no such dynamical symmetry in these shells. 

5.7 Universality of Normal and Exotic States 

The measured electromagnetic transition strengths 7? (772) between the ground states and 
first 2 + states of even-even rare earth nuclei vary with mass number A in a way quite similar 



to that of the summed orbital B (Ml) strengths measured in the same nuclei |]46 |. The 
observed correlation in these strengths suggests a microscopic relationship between modes 
that are qualitatively different in the simplest geometrical picture. Assuming that all Ml 
strength has been detected experimentally, it is difficult to explain simply within the IBM 
F-spin symmetry the observation that both the 772 and Ml strengths saturate well before 
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midshell |fl7 |. The FDSM is known to reduce to an IBM model in the limit of neglected Pauli 
correlations ||, so it is reasonable to ask whether a model such as the FDSM that incorporates 
the fermionic nature of the correlated pairs can account for the similar behavior of E2 and 
Ml strengths. 

As noted in the preceding section, the F-spin algebra of IBM-2 allows the introduction 
of an F-spin scalar (the Majorana interaction) that is insensitive to the normal symmetric 
states. Thus, It may be chosen to have a strength that pushes the mixed-symmetry states 
to their proper energies. Unlike the IBM-2, the FDSM cannot have a closed F-spin algebra. 
Thus, one is compelled to describe all states, symmetric or mixed-symmetry, in terms of a 
single Hamiltonian whose eigenstates must reflect the subtle balance between n-p quadrupole 
interactions and the n-n and p-p pairing forces. 

5.7.1 Hamiltonian and Operators 

We employ a 5-parameter (pairing and QQ) FDSM Hamiltonian: 

H = + G' 0u StS u + B' 2 ,Pl ■ Pi + B' 2v Pl ■ Pi + B 2nu Pl • P 2 V . (12) 

The operators are defined in || and the FDSM quadrupole pairing interaction is taken into 
account by renormalizing the parameters: G' 0a = Go a — G20 and B' ra = B ra — G20 (c = n,v). 
The five parameters in Eq. (O) were determined numerically using a gradient search with 



the FDUO code P7|, fitted to the experimental energies of the available 2f, 4f, 6^, 
1 + , 2 2 , and 0^ states in Nd, Sm, Gd, Dy, and Er. These calculations indicate that the 
QQ-interaction plays the crucial role in correlating the energies of the 2f and 1 + states. 
Using these parameters, we may then ask whether there is a correlation between symmetric 
and mixed-symmetry states caused by the same Q w ■ Q u strength for the electromagnetic 
transitions. 
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5.7.2 Spectrum and Multipole Strengths 



In Fig. [16| the experimental and calculated energies of the 2J~ and 1^ state energies are plotted 
versus the P factor [48], |46|, P = N p N n /(N p + N n ), where N p (N n ) are the valence proton 
(neutron) numbers, respectively. The B(E2) values and the summed B(M1) strengths 
resulting from this calculation are compared with measured quantities in Fig. |T7[ The curves 
are the empirical relations presented in [46| that summarize the approximate behavior of the 



data. The B(E2) and B(M1) strengths are reproduced quantitively by the calculations 
(with the single exception of the 164 Dy point). Thus, we find theoretical evidence for the 
approximate universal behavior of E2 and Ml strengths exhibited by the data. Furthermore, 
even the deviations from universality exhibited by the data are reproduced by the calculations 



without parameter adjustment. In Figs. |T7| c and yip, we have plotted the ratio E(A{) / E(2~l) 
as a function of P. This quantity also varies with P in a way similar to that of the E2 and 
Ml strengths, and is also quantitively reproduced by these calculations. 

Finally, we address the question of how important the variation of effective interaction 



parameters is to the success of the calculations we have described. In Figs. 17 a"— c" we 
repeat the calculations of Figs. [TTja'-c', but with a fixed set of parameters for all nuclei: 
G' 0n = -0.074 MeV, G' 0u = 0.020 MeV, B' 2w = -0.001 MeV, B' 2v = 0.047 MeV, and B 27rv = 
—0.243 MeV. The results of this calculation are seen to differ in minor details only from 
the previous calculations in which the effective interaction was permitted to have a weak A 
dependence. Thus, the quantitative reproduction of E2 and Ml strengths in the rare earth 
nuclei is an inherent feature of the FDSM with a constant effective interaction. 
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6 Higher Heritage Configurations 

The FDSM is a truncated shell model with symmetry-dictated S and D pairs. In its simplest 
implementation, one restricts the space to heritage u = 0; then all nucleons form S or D 
fermion pairs: 



4 r 



fii/2 



u ki u ki 



Or 

, r = 0, 2 for i-active 

Ofj, 

(13) 

, r = 0, 2 for k-active 



where i = ~ or k = 1, S>= and = A^, and f2i is the degeneracy for the normal-parity 
levels of a major shell. 

As we have seen, this space gives a good description for even-even nuclei in states of low 
angular momentum. However, numerical calculations within the u = space consistently 
overestimate the energies and underestimate the B(E2) values for higher angular momentum 
states, and these discrepancies increase with angular momentum. Some improvement can 
be gained in the energies by correcting perturbatively for the influence of pairing in the 
same manner as the 5*0(5) + pairing approach discussed previously. However, it has been 
realized since the inception of the FDSM that the primary reason for these discrepancies is the 
significant contribution of broken pairs to high angular momentum states [^, |50[ . Broken 
pairs in normal or abnormal parity orbitals may carry significant angular momentum, thereby 
reducing the amount carried by the S-D condensate. As was demonstrated schematically in 
Refs. |49|, 5D| , this leads to a Variable Moment of Inertia (VMI) behavior for the moment of 



inertia similar to that observed experimentally, and brings calculated and observed high-spin 
B(E2) values into quantitative agreement with data |50], [H], f|. 

Thus, it is natural to extend the numerical implementation of symmetry dictated trun- 
cation by using the FDSM to include unpaired particles. This has been discussed for a single 



unpaired particle in Refs. [^2|, |53, B4|, || . We shall not discuss that further here, but instead 
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conclude by summarizing recent work that incorporates broken pairs in such calculations. 
We have developed a computer code SU3su2 that includes explicit broken pairs in an SU(3) 
limit in order to examine the yrast states in the rare earth and actinide regions [[5^, |5(J. In 
this extension, the model basis is constructed by coupling the SU(3) basis and one broken 
pair in the normal parity levels or the unique parity level 

(S, D) N ^ <g> A'(i) or (S, D) N ^ <g> A'(j ). (14) 

where S and D are the usual symmetry-dictated coherent fermion pairs with coupled angular 
momentum and 2, respectively, A'(i) designates broken normal-parity pairs, and A' (jo) 
designates broken pairs in the unique parity level. The corresponding creation operators are 

(KI)L' 



A't(i) = &L&L ' > (is) 



hi 



^ t (jo)=EA/ 2 No^ t o]l; • (is) 



ki 

The basis for the SU(3) core is 

Sp(6) D SU(3) D SO(3) 
Ni, u = (A, ji) kL 

where k is an additional quantum number to distinguish orthogonal states having the same 
(A,/i) and J. For an SU(3) core and one broken pair, the following coupling schemes are 
possible: 

\N lU = 0(Ai, fix) ® (2, 0); (A 2 , y^KIJMj), 

\N 1 u = 0(\ 1 ,n 1 )KL®K i ;KIJMj), (18) 
\N x u = 0(Ai, h x )kL ® I ; JMj). 

and the total Hamiltonian is 

H = H SD + H A > + H mix , (19) 
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where ifsD is the S£7(3)-plus-pairing Hamiltonian, Ha> corresponds to the broken pair, and 
-f^mix is the interaction between the (S, D) core and the broken pair. The last term H m i x 
leads to the mixing of heritage. Explicitly, the Hamiltonian can be expressed as 

H = aK-K + jI-I-5I-K + j I -I -5 Q Io-K 

- KP 2 {k) ■ P 2 (i) + k'(D^D) 2 • P 2 (i) - K P 2 (k) ■ P 2 (i ) + n' (D^D) 2 ■ P 2 (i ) 

+ AB , (20) 

where AB is the energy required to break one pair in either normal parity or abnormal parity 
levels. 



In Fig. we present the results of a u = 2 calculation of the energy levels up to J = 20 
for the even-even isotopes 160-166 Er |56||. The inclusion of a single broken pair in the ST/3 
symmetry limit leads to a quantitative description of the spectrum. We do not show results 
for B(E2) values from these calculations because that work is still in progress. Preliminary 
indicatons are that the inclusion of a single broken pair leads to a substantial increase in 
the B(E2) strengths at higher spins, but the calculated values are still low in the angular 
momentum 10-20 region, suggesting that configurations having two broken pairs may be 
important for a quantitative description of B(E2) values near angular momentum 20 for 
rare earth nuclei. 

7 Conclusions 

The Fermion Dynamical Symmetry Model (FDSM) provides a systematic method for trun- 
cating the spherical shell model. In this scheme, a valence space is selected using energy 
considerations and principles of dynamical symmetry are then used to radically truncate the 
valence space. We term this a symmetry-dictated truncation. The resulting truncated space 
permits systematic shell model calculations to be implemented for all heavy nuclei. Since 
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the space has been severely truncated, the corresponding interactions are highly effective 
with respect to the original shell model. Thus, the first step in systematic calculations with 
this truncation scheme is to determine the appropriate FDSM effective interaction for each 
valence space. Although it is of considerable interest in the longer term to relate such effec- 
tive interactions to standard shell model ones, the most practical initial way to determine 
the required interaction is to construct it phenomenologically by viewing its matrix elements 
as parameters to be constrained by a carefully chosen data set. We have presented examples 
of systematic FDSM calculations that have been used to determine an effective interaction 
appropriate for configurations in heavy nuclei with no broken pairs. This interaction appears 
to be simple and to have a rather weak dependence on particle number within major-shell 
valence spaces. Calculations using this interaction reproduce low-lying spectra, moments, 
and transition rates for broad ranges of nuclei that exhibit varied collective behavior: axial 
rotors, anharmonic vibrators, and gamma-soft rotors. Finally, we have presented an initial 
extension of this approach to include broken pairs in the configuration space. These results 
constitute a practical demonstration that systematic shell model calculations are now fea- 
sible for very heavy nuclei far removed from closed shells. The primary remaining task is 
to systematize and refine the interaction over all mass numbers through a detailed series of 
numerical calculations. 
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(Shell Model Truncations) 
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Example: Energy-Dictated Truncation 
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Example: Symmetry-Dictated Truncation 
H(SU^U,)~af + bJl 
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Spontaneous Symmetry Breaking 

Phase Transitions 
V J 



Figure 1: An illustration of energy-dictated and symmetry-dictated trun 
cations for a simple symmetry. 




Energy-Dictated Truncation 




Symmetry-Dictated Truncation 
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Figure 2: Energy-dictated and Symmetry-dictated truncations for a spher 
ical shell model. 
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Figure 3: The FDSM coupling scheme for normal and superdeformation. 



Dynamical Symmetry: 

G, => G 2 3 G 3 -- G„ H sym = f(C G ) 

Sp 6 zd SU 3 zd S0 3 H = H Q + aC Sp(6) + bC su(3) + cC SO(3) 



Dynamical Symmetry Breaking: 

/ ZJ \ = /ff\_/ff \ 




= F -'hi \ 

H=H sym + H^ h < \ M sym/ 



(Empirical) 



H sb = ^micro - H sym (Microscopic) 



Figure 4: Dynamical symmetry and the breaking of dynamical symmetry. 
The Casimir invariants of the group G are denoted by Cq. 
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QUESTIONS 

1. How large are the symmetry 
breaking terms H sb ? 

2. What is the relevant scale 
against which to measure the 
symmetry breaking terms? 

3. How much of the deviation of 
the empirical energies from the 
energies assumed in the 
symmetry limits actually breaks 
SU(3)? 



Figure 5: The difference between the symmetry- limit and realistic single- 
particle spectrum. 




Figure 6: Range of upper limits for SU3 symmetry breaking by single- 
particle energies in the 82-126 shell. 
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Figure 7: Mass shell correction for heavy and superheavy elements using 
the dynamical symmetry methods of Ref. |36[| , 
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Figure 8: Comparison between calculated levels using eq. (||), and experi- 
mental energy levels for the even-even 120_126 Xe isotopes. 
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Figure 9: Comparison between calculated levels using eq. (|8|) and experi- 
mental energy levels for the even-even 128_132 Xe isotopes. 
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Figure 10: Comparison between calculated levels using eq. @ and exper- 
imental energy levels for the even-odd 127-129 Xe isotopes. The even-odd 
nuclei are constructed by coupling the neighboring even-even core to a va- 
lence neutron. 
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Figure 11: Comparison between calculated levels using eq. @ and exper- 
imental energy levels for the even-odd 131_133 Xe isotopes. The even-odd 
nuclei are constructed by coupling the neighboring even-even core to a va- 
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Figure 12: The systematic shift of the ground band as a function of mass 
number for Xe isotopes. 
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Figure 15: Spectra for positive-parity states in 190 194 > 198 pt. Experimental 
levels are taken from [58, 5£, |6l| ]. 
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Figure 16: Experimental and calculated energy levels for the first 2 + 
(lower points) and 1 + (upper points) states in selected even-even rare earth 
nuclei where orbital Ml strengths have been measured. The symbols denote 



isotopes with the same meanings as in Fig. 17 



54 



M. W. Guidry, D. H. Feng, X.-W. Pan, and C.-L. Wu 



03 



C\2 

m 



15 



10 











■ 1 


1 1 1 1 1 1 1 1 


' 1 


o 


168 Er 


Expt. 




O 


i 8 a Er 


Theo. 


- o 


160-164 Dy 




10 


- O 


160-164 Dy 






154-160 Gd 




-k 


154-160 Gd 




A 


144,148-154g m 






A 


144,148-1542 m 




□ 


142,148-150 Nd 




5 


□ 


142,146-180^^ 








* (a) 








f (a') 


. -A- 






n 






. i ...... . 







5 


c 







5 



PQ 



+ 



+ 




15 
10 

- 5 


6 



- 2 



1 1 

; o 


i | i | i | i 


i ' 1 

1 11C \J . 


- o 


160-164 Dy 






154-160 Gd 




A 


144,148- 154c™ 




□ 


143,146— 150^ ^ 












□ 


API / i i \ ' 

y (a") 

.. 1 ...... . 







5 


1 i 

■ 


-. — | — i — | — i — | — i- 

168 Er 


' ■ i 

Theo. . 


o 


160-164-Qy 




- -k 


154-160 Gd 




A 


144,148-1542 m 




- □ 


142,146- 150j^ ^ 














A X A X 




A 






□ 


7 ( b "> 




j i i i i i ■_ 





— 1 — i — 1 — 


i — i — i — i — i — > — i — i — i — ' — i — i — i — i — i — i — 


u 


1 1 1 


i — i — i — ' — i — i — i — ' — | — i — i — i — i — i — i — i — 


u 


1 1 ' 


1 — 1 — 1 — ' — 1 — ' — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — ' — 




Expt. . 






Theo. . 






Theo. . 






4 






4 




□ 

a a ^dfcfo ■ 


a — 


n/4 


2 






2 


a 


^U^U^. 




(c) 






(C) 






(c"> 









. 1 . 


1 





. i . 


1 





5 

P=N p N n /(N p +N I1 ) 





5 

P=N p N n /(N p +N n ) 





5 

P=N p N n /(N p +N„) 



Figure 17: Comparison of experimental and theoretical B(E2) and 
summed B{M1) values for rare-earth nuclei as a function of P. The ex- 
perimental results in the left column are taken from the compilation given 
m Q and p§. For comparison, curves from the same empirical relation 
employed there, B(E2, Ml) = a\ + «2/[l + exp((c — P)/d)} are also plot- 
ted. The symbols in (c)-(c") have the same meaning as in the E2 and Ml 
cases. The theoretical results in the second column of figures correspond 
to best adjustment of effective interaction parameters to reproduce spectra. 
The theoretical results in the third column of figures correspond to constant 
values of the effective interaction parameters for all nuclei examined. 
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Figure 18: Experimental and theoretical energy levels in some even-even Er 



isotopes calculated with a single broken pair using the code SU3su2 J55], 56]. 
The spectrum of 160 Er is to the left and that of 166 Er to the right. 
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Table 1: Experimental and Theoretical Mass Excesses in MeV for the Heaviest Elements. 



256 104 


258 105 


260 106 


262 107 


264 108 


266 109 


269 U0 


272 ln 


Reference 


94.25(03) 


101.84(34) 


106.60(04) 


114.68(38) 


119.82(30) 


128.39(35) 


134.80(32) 


141.70(37) 


Exp |Jsa 


J35| 


94.15 


101.98 


106.8 


114.82 


120.02 


128.16 


133.31 


140.49 


FDSM 




95.90 


103.41 


106.87 


116.10 


121.28 


129.44 




144.83 


Myers 1 


34F 


94.84 


102.22 


105.81 


115.00 


120.40 


128.43 




144.04 


Groote 





95.6 


102.6 


106.8 


114.9 


120.4 


127.8 




142.6 


Scegcr J 




94.37 


101.64 


105.68 


114.78 


120.27 


128.44 




143.82 


Liran [p 




95.77 


103.01 


108.13 


115.50 










Moller 


3q| 


95.69 


103.11 


108.12 


115.71 


121.09 


129.04 


135.51 


143.44 


Moller 


36J 


93.78 


101.00 


105.81 


113.18 


118.34 


126.06 


132.39 


140.18 


Moller 




93.38 


100.97 


105.73 


113.45 


118.73 


126.65 


133.08 


140.93 


Moller 


3q| 


94.52 




107.04 




120.47 








Patyk [ 




94.38 




106.93 




120.47 




135.46 




Cwiok 


70) 



Solution of the Nuclear Shell Model by Symmetry-Dictated Truncation 



Table 2: Parameters for even Xe isotopes. 



Nuclei 


ffe(keV) 


^(keV) 


^(keV) 


120 Xe 


-60 


53 


11.9 


122 Xe 


-64 


59 


11.9 


124 Xe 


-68.8 


64 


11.9 


126 Xe 


-73.3 


71 


11.9 


128 Xe 


-78.2 


79 


11.9 


130 Xc 


-100.9 


100 


11.9 


132 Xe 


-106.5 


122 


11.9 



Table 3: Parameters for odd Xe isotopes. 



Nuclei 


<? 6 (keV) 


^(keV) 


5 5 (keV) 


Sj(keV) 


127 Xe 


-73.3 


71.0 


-38.0 


25.0 


129 Xe 


-78.2 


79.0 


-18.0 


35.3 


131 Xe 


-100.9 


100.0 


30.0 


35.3 


133 Xe 


-106.5 


122.0 


50.5 


35.3 


135 Xe 




142.1 


70.6 


35.3 



Table 4: Parameters for odd Ba isotopes. 
Nuclei 5 |(keV) 55 (keV) <?j(keV) 



131 Ba 


72.5 


-20 


35.3 


133 Ba 


105 


-15 


35.3 


135 Ba 


90 


50 


35.3 


137 Ba 


72 


70 


35.3 
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Table 5: Relative B(E2) values for the even Xe isotopes. 







120 


Xe 


124 


Xe 


126 


Xe 


128 


Xe 


130 


Xe 


J* 




Exp. 


Theo. 


Exp. 


Theo. 


Exp. 


Theo. 


Exp. 


Theo. 


Exp. 


Theo. 


2+ 


-2+ 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 




-of 


5.6 


5.6 


3.9 


3.9 


1.5 


1.4 


1.2 


1.2 


0.6 


0.6 




^2+ 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 




-4+ 


50 


40 


46 


40 


34 


40 


37 


40 


25 


40 




-2+ 


2.7 


7.1 


1.6 


4.9 


2.0 


1.85 


1.0 


1.5 


1.4 


0.72 


4 


^2+ 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 




-4+ 


62 


91 


91 


91 


76 


91 


133 


91 


107 


91 




-2+ 




7.11 


0.4 


4.91 


1.0 


1.83 


1.7 


1.49 


3.2 


0.97 


0+ 


^2+ 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 




-2+ 




7.11 


1 


4.91 


7.7 


1.83 


14 


1.49 


26 


0.97 
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Table 6: Transition probablities in 129 Xe (JVi = 5) and 131 Xe (iV 1= 4), Q 



20. 



129 Xe 



131 Xe 



B{E2)(e 2 b 



B(E2){e 2 b 



Ji J f FDSM Exp. Ref. Q 



Ji -» J/ FDSM Exp. Ref. [[T| 



3 + 
2i 


_^ 1 + 
2i 


0.036 




0.007 


3 + 
2 2 


. 3 + 
2i 


0.0186 


<0.0005 


0.013 


3 + 
2 2 


_^ 1 + 
2i 


0.084 


0.12 


0.12 


5 + 
2i 


^ 3 + 
2i 


0.011 


0.22 


0.10 


5 + 
2i 


_^ 1 + 
2i 


0.070 


0.077 


0.039 


1 + 
2 2 


, 3 + 
2 2 


0.028 






1 + 
2 2 


. 3 + 
2i 


0.14 


0.044 


0.12 


1 + 
2 2 


_^ 1 + 
2i 


0.0000 






5 + 
2 2 


_^ 1 + 
2i 


0.004 


0.057 


0.071 


3 + 
2 3 


_^ 1 + 
2i 


0.056 


0.0032 


0.0056 


3 + 
2 4 


_^ 1 + 
2i 


0.0133 


0.0030 


0.0004 



1 + 

2i 

5 + 
2i 

r ± + 
2i 

3 + 
2 2 

1 + 
2 2 
1 + 
2 2 

1 + 
2i 

Z+ 

2i 

3 + 

2 3 

5 + 
2 2 

2 2 

^2 + 
2 2 

1 + 

2, 



3 + 
2i 

1 + 

2 a 

3 + 
2 a 

3 + 
2i 

1 + 
2 a 

3 + 
2 a 

5 + 
2 a 

3 + 
2i 

3 + 
2 a 

3 + 
2 2 

1 + 

2 a 

3 + 
2 a 

3 + 
2i 



0.0953 

0.075 

0.004 

0.053 

0.0000 

0.124 

0.028 

0.043 

0.025 

0.004 

0.071 

0.014 

0.124 



0.0039 0.0012 
0.030 0.016 



0.10 
0.057 

0.048 
0.005 
0.081 
0.027 
<0.031 
0.068 
0.013 
0.005 



0.10 
0.058 

0.115 

0.0013 

0.082 

0.017 

0.0011 

0.056 

0.043 

0.026 



M. W. Guidry, D. H. Feng, X.-W. Pan, and C.-L. Wu 



Table 7: Parameters for the Pt calculations. 



Mass= 


190 


192 


194 


196 


198 


N lu 


6 


6 


5 


4 


3 




-386 


-386 


-386 


-386 


-386 


B2tt 


48 


48 


48 


48 


48 


Gq v 


-49 


-49 


-49 


-49 


-49 


Gon 


-10 


-3 


-5 


-18 


-50 


B2u 


36 


53 


74 


97 


100 


e n 


0.17 


0.17 


0.15 


0.16 


0.18 




0.19 


0.19 


0.15 


0.14 


0.19 




g n = 0.62, 




9u = 








13 x 10- 3 fm, 


(3ou = -1.2 x 10 


- 3 fm 
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Table 8: Comparison of B(E2) values for 196 Pt. 



Ji 




£(£2)^ p 


B(£2)gj, 


50(6) limits 


IBM-2( £ ) 


IBM-l(g) 


FDSM 


It 




0.288(14) 


0.274(1) 


0.288 


0.289 


0.288 


0.195 


4 




0.403(32) 


0.410(6) 


0.378 


0.395 


0.393 


0.248 






0.421(116) 


0.450(28) 


0.384 


0.409 


0.423 


0.215 


8f 


-6J- 


0.577(58) 




0.341 


0.325 


0.416 


0.139 


2+ 

z 2 


-2+ 


0.350(31) 


0.370(5) 


0.378 


0.40 


0.303 


0.262 


2+ 


->o+ 


< 2.0 x 10~ 6 W 







0.001 


0.004 


0.0001 


0+ 


^2+ 


0.142(77) 


0.1(1) 


0.385 


0.466 


0.375 


0.268 


0+ 


-2+ 


0.033(7) W 


0.028(5) 





0.026 


0.007 


0.021 


4 


-4+ 


0.193(97) W 


0.084(14) 


0.183 


0.206 


0.171 


0.121 


4 


^2+ 


0.177(35) W 


0.18(2) 


0.201 


0.206 


0.199 


0.112 


4 


-2+ 


0.0030(10) W 


0.001(2) 





0.006 


0.004 


0.001 


6+ 


-6+ 


0.085(121) W 




0.108 


0.12 


0.11 


0.075 


6+ 


-4+ 


0.350(102) W 




0.232 


0.201 


0.25 


0.056 


6+ 


-4+ 


0.0037(16) W 







0.017 


0.001 


0.031 


n 


-2+ 


0.0009(15) W 







0.004 


7.2 x 10~ 6 


0.004 


o 3 + 


-2+ 


< 0.034 ( c ) 







0.067 


0.003 


0.010 



(°) B(E2) cxp are taken from H. H. Bolotin et al, Nucl. Phys. A370 (1981)146 

W Data are from A. Mauthofer et al, Z. Phys. A336 (1990) 263. 

W Data are from H. G. Borner et al, Phys. Rev. C42 (1990)R2271. 

(<*) B(E2) cxp are taken from C. S. Lim et al, Nucl. Phys., A548 (1992) 308. 

(e 1 Calc. from Bijker et al, Nucl. Phys. A344 (1980) 207. 

V) For 0(6) limit, e = 0.155(e6). 
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Table 9: Comparison of experimental and theoretical g-factors in 190 198 Pt. 



o+ 4+ o+ 

z l *1 z 2 



isotope 


exp>) 


IBM^ 


FDSM 


exp>) 


IBM FDSM 


cxp>) 


IBM FDSM 


190 






0.137 




0.092 




0.121 


192 


0.318(17) 


0.23 


0.182 


0.281(30) 


0.126 


0.278(46) 


0.139 


194 


0.295(10) 


0.27 


0.239 


0.279(31) 


0.171 


0.281(55) 


0.178 


196 


0.295(10) 


0.32 


0.268 


0.345(40) 


0.213 


0.271(45)^ 


0.207 


198 


0.293(34)( b ) 


0.40 


0.26 


0.307(54)W 


0.170 


0.307(54)( b ) 


0.213 



Data are taken from F. Brandolini et al., Nucl. Phys. A536 (1992) 366. 

W Data are taken from A. E. Stuchbery, G. J. Lampard and H. H. Bolotin, 
Nucl. Phys. A365 (1981) 317; A528 (1991) 447. 

( c ) IBA-2 calc. from M. Sambataro, et al, Nucl. Phys. A423 (1984) 333. 



Table 10: Experimental and theoretical changes of charge radii for 190 198 Pt isotopes. 







<5 < 


r 2 > (10~ 3 


fm) 


A < r 2 > (10~ 3 fm) 




190 


192 


194 


196 198 


192-190 


194-192 196-194 


198-196 




2+ 2+ 
z l z 2 


2+ 2+ 
z l z 2 


2+ 2+ 
z l z 2 


2 + 2 + 2 + 2 + 
z l z 2 z l z 2 




FDSM 


4.48 3.22 


3.84 4.76 


3.16 4.93 


2.66 5.54 2.14 3.68 


74 


75 71 


69 


IBM 






3.62 5.08 


2.55 








Exp» 






3.45 3.57 


4.49 


65 


67 74 


80 



(a) Isomer shifts arc taken from R. Engfer et al, Nucl. Data Tables, 14 (1974) 509; 
Isotope shifts are taken from Th. Hilberath et al, Z. Phys., A342 (1992) 1. 



